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ABSTRACT 

Recent observations of Cepheids in the Virgo cluster have bolstered the evidence that supports a 
Hubble constant in 70— 90 km sec -1 Mpc -1 range. This evidence, by and large, probes the expansion 
of the Universe within 100 Mpc. We investigate the possibility that the expansion rate within this 
region is systematically higher than the true expansion rate due to the presence of a local, large 
underdense region or void. We begin by calculating the expected deviations between the locally 
measured Hubble constant and the true Hubble constant for a variety of models. The calculations 
are done using linear perturbation theory and are compared with results from N-body simulations 
wherever possible. We also discuss the expected correlations between these deviations and mass 
fluctuation for the sample volume. We find that the fluctuations are small for the standard cold 
dark matter as well as mixed dark matter models but can be substantial in a number of interesting 
and viable nonstandard scenarios. However, deviations in the Hubble flow for a region of radius 
200 Mpc are small for virtually all reasonable models. Therefore, methods based on supernovae or 
the Sunyaev-Zel'dovich effect, which can probe 200 Mpc scales, will be essential in determining the 
true Hubble constant. 

We discuss, in detail, the fluctuations induced in the cosmic background radiation by voids 
at the last scattering surface. In addition, we discuss the dipole and quadrupole fluctuations one 
would expect if the void enclosing us is aspherical or if we lie off-center. 
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1 Introduction 



It is straightforward to measure the Hubble constant in a perfectly homogeneous and isotropic 
Universe: one simply divides the relative velocity between any two points by their separation. The 
situation is of course far more difficult in our inhomogeneous and anisotropic Universe. For most 
cosmologists, the working assumption is that the Universe becomes homogeneous on very large 
scales. The goal then is to find the recessional velocity and distance to an object sufficiently far 
away so that the Hubble flow dominates over peculiar velocities. In this work, we address the 
possibility that determinations of the Hubble constant have not yet reached such distances. 

Knowledge of the true, global value of the Hubble constant, Hq, is central to any cosmological 
model. Most measurements to date place Hq between 65 — 95 km sec -1 Mpc" 1 , values that are 
uncomfortably high given estimates of the age of the Universe based on globular cluster dating, and 
the theoretical prejudice that we live in a spatially flat, matter-dominated universe. In particular, 
if we accept that the Universe is older than 12 Gyr, then a Hubble constant h > 0.55 (Hq = 
100 h km sec -1 Mpc -1 ) is incompatible with all Q m = 1 models (J) m is the density of matter in units 
of the critical density) and h > 0.80 is incompatible with all zero cosmological constant (A = 0) 
models. In addition, a high value for the Hubble constant is, in general, trouble for structure 
formation scenarios. This is especially true for cold dark matter (CDM) where the favoured value 
is h ~ 0.3 — 0.4. Finally, a low value for Hq helps to close the gap between the prediction for the 
baryon density from big bang nucleosynthesis and recent determinations of the dark-to-baryonic 
matter ratio in rich clusters (Bartlett et al. 1994). 

All resolutions to the "age problem" come at a cost. A true age for the Universe significantly 
less than 12 Gyr would require radical rethinking of stellar evolution (Cage & Walker 1992; Sandage 
1993; Chaboyer 1994; Shi 1995). On the other hand, we can accommodate a universe older than 
12 Gyr by requiring the Universe to be open, or by invoking a non-zero cosmological constant. For 
fixed h, an open Universe corresponds to an older one; but this is problematic especially if one 
accepts inflation which predicts O = 1 in nearly all incarnations (see however Bucher et al. 1995, 
Linde 1995). Similarly, a universe dominated by a cosmological constant (or any energy density 
with negative pressure) will be older than the corresponding f2 m = 1 universe. The cost here is the 
introduction of a new energy density term in the Friedmann equations whose origin is completely 
mysterious. 

In many respects, a true Hubble constant of 50kmsec _1 Mpc -1 or less represents the simplest 
solution to the age problem. However, the majority of observations suggest that such values are 
ruled out. These include recent measurements by the Canada-France-Hawaii Telescope (CFHT) 
(Pierce et al. 1994) and the Hubble Space Telescope (HST) (Freedman et al. 1994) whose published 
values of h = 0.87±0.07 and h = 0.80±0.17, respectively, exclude a spatially flat, matter dominated 
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universe. Their Cepheid results are also consistent with other determinations based on techniques 
such as Tully-Fisher, planetary nebulae, and surface brightness fluctuations (see van den Bergh 1992 
for a review). In principle, Cepheids should provide the most accurate method for determining the 
distance to Virgo. At present, there is the potential for large systematic errors (see, in particular, the 
paper by Freedman et al. 1994) due primarily to uncertainties in the position of the galaxy studied 
within Virgo. Measurements of Cepheids in other galaxies within the cluster will significantly 
reduce this uncertainty and provide a distance determination to Virgo of better than 10%. 

Virgo is only 17Mpc away, far too close for a reliable measurement of Hq. But it is only one 
additional rung of the distance ladder to clusters such as Coma which, at a distance ~ 100 Mpc, are 
assumed to be far enough away that their recessional velocities are dominated by the true Hubble 
flow. 

But is this last assumption valid? Or is it possible that the local Hubble flow, as measured out 
to ~ 100 Mpc, overestimates the true Hubble constant by 20-30%? The discovery of large voids 
and sheet-like structures in redshift surveys (Kirshner et al. 1981; de Lapparent et al. 1986; Geller 
& Huchra 1989; da Costa 1991) and large scale bulk flows (Lauer & Postman 1994) on ~ 100 Mpc 
scales should give us pause as variations in the Hubble flow are likely to occur on comparable scales. 

The question of deviations in the Hubble flow was addressed by Turner, Cen, and Ostriker 
(1992, hereafter TCO) in the context of a simple numerical experiment. A model universe (with a 
specified value for the global Hubble constant) is simulated using standard N-body methods and 
it is assumed that "observers" in this universe can accurately measure the distances to and radial 
velocities of all galaxies within a set distance R. Each observer constructs a Hubble diagram and 
reads off a value for the local Hubble constant . In this way, the effects of local deviations in the 
Hubble flow are separated out from observer-dependent effects such as systematic errors in distance 
determinations and incomplete sampling. The distribution of Hi, measured by all observers in the 
simulation is studied and compared with other statistics of the density and velocity fields such as 
rms mass fluctuations and bulk flows. In the end, one can make probabilistic statements about the 
likelihood a given observer has for measuring a local Hubble constant within the range allowed by 
observations. 

An alternate approach to this global statistical method is to explore specific and simple models 
for our local region of the Universe. Suppose, for example, we live in an underdense region or 
void. Any determination of the Hubble constant based on objects within this region will reflect 
the local expansion rate rather than the global one. The simplest model of this type has us at 
the center of a spherically symmetric section of an open spacetime embedded in a spatially flat 
Friedmann- Robertson- Walker universe (Suto et al. 1994; Wu et al. 1995; Moffatt & Tatarski 1994). 

In Section 2 we take a second look at the statistical approach pioneered by TCO. Our interest is 
in studying a variety of models and so we use linear perturbation theory where probabilities can be 
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calculated directly from the linear power spectrum P(k). Our first step is to compare linear theory 
results with those of TCO. As expected, agreement is good on large scales where perturbations 
are small. On small scales, linear theory tends to underestimate the rms fluctuations in Hi. In 
addition, TCO find that the distribution of measured Hi has a distinctly non-Gaussian shape. We 
next use linear theory to study correlations between Hi and other observables. In particular, we 
derive the statistical relationship between Hi and the mass excess for the sample volume. 

TCO consider the standard CDM scenario circa 1992 which is now known to be inconsistent 
with results from the COBE DMR experiment (Smoot et al. 1992) and the APM galaxy survey 
(Maddox et al. 1990). Here we use linear theory to calculate the rms fluctuations in Hi for a 
variety of observationally viable models. As one might expect, the size of deviations in the Hubble 
constant on 100 Mpc scales depends sensitively on the shape of the primordial power spectrum but 
is rather insensitive to the type of dark matter (i.e., hot, cold, or mixed) in the model. 

Section 3 deals with the more direct approach of modeling our local region of the Universe. The 
simplest scenario is to assume that we are near the center of an underdense region that extends 
well past the Coma cluster. The mean density in this region must be fairly low (f2 vo id ^ 0.5) in 
order that there be a significant difference between the true and measured values for the Hubble 
constant. This may not be so farfetched: there are hints from the number counts of galaxies that 
we live in a very large underdense region (Metcalfe et al. 1992) though interpretation of this data 
is complicated by evolutionary effects (see, for example, Loveday et al. 1992). 

A large underdense region today implies rather significant perturbations in the energy density 
at the surface of last scattering. These perturbations lead to fluctuations in the cosmic back- 
ground radiation (CBR) through the Sachs-Wolfe effect. We can therefore use CBR anisotropy 
measurements to limit the number of voids that reside in our Hubble volume. In addition, large 
angular-scale perturbations will be induced if our local underdense region is aspherical, or if we lie 
off-center. We can therefore place important constraints on this scenario by using the dipole and 
quadrupole anisotropy measurements. 

Section 4 summarizes our results. Our general conclusion is that large fluctuations in Hi for 
R ~ 100 Mpc can occur, but only in nonstandard models. For example, a CDM model with a tilted 
primordial power spectrum P{k) oc k n ;n > 1.5 will have 2a fluctuations in H in excess of 30%. 
However, even in these nonstandard models, fluctuations in H fall rapidly with increasing scale: 
Determinations of H based on objects ~ 200 Mpc away should reflect the global value to within 
10% unless we live in an extremely rare, very large-scale underdense region. 

There is, at present, considerable discrepancy among measurements for objects > 200 Mpc 
away. The "expanding photosphere method" has been used to determine distances to 18 type II 
supernovae as far away as 180 Mpc (Schmidt et al. 1994) and references therein). An analysis of 
this data finds h = 0.73 ± 0.06 (statistical) ± 0.07(systematic) km sec -1 Mpc -1 with no evidence for 
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spatial deviations in the Hubble flow. In addition Lauer & Postman (1994) find that the value of 
H does not vary by more than 7% between 30 and 150/i _1 Mpc. On the other hand, measurements 
for very distant (R > 400h _1 Mpc) objects based on the Sunyaev-Zel'dovich effect do seem to yield 
systematically lower values for h (Yamashita 1994; Birkinshaw & Hughes 1994) suggesting that we 
may in fact live in a low density, high h region of the Universe (Suto et al. 1994). The supernovae 
results and Lauer & Postman data seem to indicate that this void would be much larger than 
anything expected, even for fairly nonstandard models designed to give more power on large scales. 

2 Hubble Constant Statistics 
2.1 Linear Theory Analysis 

In linear perturbation theory (Peebles 1993), the relationship between the peculiar velocity field 
and the density contrast <5(x) = (p(x) — pb)/pb is given by 

where /o(x), pi, and Hq are the density field, background density, and global Hubble constant today. 
The peculiar Hubble flow, 5H(x) = Hl(x.) — Hq is 

6H{x) = [ d 3 yv(y) ■ JL^-H/( X - y) (2) 

J |x-y| 

where Hl(x.) is the local value of the Hubble constant as measured by an observer at position x and 
W(x— y) is the window function for the observations. 5H is essentially a breathing mode (expansion 
or contraction) of the peculiar velocity field, corresponding to the trace of the shear tensor (see, 
e.g. Kaiser 1988). For a perfectly volume limited sample out to a distance R, W is a step function, 
W(x — y) = 9(R — |x — y|)V[^ where Vw = 4-7ri? 3 /3. In a Friedmann-Robertson- Walker universe, 
5H -> for R -> oo. 

Equations (|l|) and (||) can be written in terms of the Fourier transform of the density contrast: 
()// ; ' d '' k 6(k)Z(kR)e^ (3) 



#o J (2?r) 3 / 2 
where 5(k) = (2vr)~ 3/2 / d 3 ^(x)e ik x , 

Z(x) = 3 Sln(x) = Si(x) , (4) 

and Si(x) = Jq dxsin(x)/x. The rms fluctuations in the local Hubble constant 5h = {{5H / Hq) 2 ) 1 / 2 
((...) denotes spatial average) is given by 

5% = —s / k 2 dkP(k)Z 2 (kR) (5) 
2tt z J 

where P(k) = \5(k)\ 2 is the power spectrum. If the primordial perturbations are Gaussian, the 
distribution of 5h will be Gaussian, at least in the linear regime. 
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R (Moc) 


5h (TCO) 


5h (XT) 


m (LT) 




10 


1.08 


0.42 


-0.64 


1.07 


20 


0.45 


0.23 


-0.64 


0.79 


40 


0.18 


0.11 


-0.65 


0.55 


60 


0.10 


0.071 


-0.67 


0.43 


80 


0.062 


0.049 


-0.67 


0.35 


120 


0.029 


0.027 


-0.67 


0.26 



Table 1: TCO and Linear theory results for 5jj 

2.2 TCO and Linear Theory 

TCO use numerical simulations of cosmological models to generate a set of hypothetical observers. 
These observers can accurately measure the distances to and velocities of the galaxies within a 
pre-set distance. The local value of the Hubble constant, as measured by the fcth observer, is given 
by 

1 x-^ r • Vj 

where rj and v; are the position and velocity vectors of the ith galaxy, r = r; — r^, and the sum is 
over all ./V galaxies in the prescribed volume. 

Eq @ is essentially the discrete version of Eq (||) with the important distinction that Eq (||) 
involves a volume average whereas Eq (||) involves a sum over galaxies. Since galaxies tend to reside 
in high density regions (this will depend on biasing, i.e., galaxy formation) where the expansion 
rate is slower, the results of TCO will tend, on average, to be skewed towards negative values of 
5H = (Hh — Hq) /Hq. This is indeed found to be the case. 

Column 1 of Table 1 summarizes the results for 5h from TCO. The results derived from linear 
theory for their model (baryon density, in units of the critical density, £Ib = 0.05; h = 0.5; f2 m + 
£Ib = 1; rrns mass fluctuation on 8/i _1 Mpc, as = 0.67) are shown in Column 2. As expected, 
the agreement is best on large scales where perturbations are small. On small scales, linear theory 
tends to underestimate 5h- In any case, it is apparent that for this particular model, the local 
value of the Hubble constant, as measured by observers with complete coverage out to 100 Mpc, is 
expected to be with 5% or so of the true or global value. 

2.3 Correlations between 6 H and other observables 

A volume limited measurement of H cannot separate the peculiar Hubble flow within that volume 
from the true Hubble flow. There are, however, other quantities accessible to observers which char- 
acterize the peculiar velocity field and matter distribution. Correlations between such quantities 
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and the peculiar Hubble flow, if they exist, would provide a means for correcting our measurements 
of H L . 

An obvious quantity to consider is mass fluctuation within the sample volume, 5M/M. For an 
observer at the center of a spherically symmetric region with constant density contrast, 5H/Hq = 
— T^SM/M so long as S <C 1. If this were the case, knowledge of 5M/M would be sufficient 
to correct local measurements of H. In general, however, there is only a statistical relationship 
between 5H/Hq and 5M/M. TCO explore this relationship by plotting 5M/M versus 5H/Hq for 
the observers in their simulation (their Figure 6a). A linear fit is obtained by minimizing the total 
orthogonal distance of all points on this figure to the straight line 
5H 5M 

w = m w + b > (7) 

that is, by minimizing the function 

Their result, m = —0.60, b = —0.008, corresponds to a slope that is nearly twice what one would 
naively expect from linear theory, and they attribute this difference to nonlinear effects. While 
they find substantial scatter in the plot, they conclude that 5M offers the best hope for correcting 
local measurements of the Hubble constant. 

Interestingly enough, the results of TCO are similar to those found in linear theory. Consider 
the linear theory expression for the fractional mass excess in a spherical top-hat region: 
SM f d 3 k „, ik . x 



M J (2tt) 3 / 2 

where 



<5(k) W{kR) e lk ' x (9) 



W(x) = 3 Sin ^-" COS ^ . (10) 
The rms mass fluctuation is given by 



S 2 M 



^2 J k 2 dkP(k)W 2 {kR). (11) 



In the limit kR — ► 0, Z — > —1/3 and W — > 1 which would seem to support the naive expectation. 
However, the window functions in Eq (|5|) and (|ll|) differ not only in amplitude but in shape with 
Z having more support at larger k. A linear fit, analogous to the one calculated in TCO, yields 



m = V 2 /SMSH\ " ~ W 

\ M H / 

where 

(WW) = k2 dkP(k)Z{kR)W(kR) (13) 
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The rms fluctuation from the best fit curve cjhm is given by 



a HM — 



6% + m 2 5\ 



:2 _i m lmm. 

M \ M H 



) 



(14) 



(1 + m 2 ) 



Our results, shown in Table 1, illustrate that those of TCO can be largely explained by linear 
theory. 

In principle, one can use knowledge of 5M/M in our region of the Universe to correct local 
measurements of H. With this in mind, we choose to determine the correlation between 5H/Hq 
and 5M/M by a slightly different procedure. Specifically, we fit minimize the distance of 5H/Hq 
to the curve given by Eq (0), i.e., we minimize 



A standard least squares analysis gives m = yjf^jf/ /$m and <r 2 HM = 5 2 H — 2m (Jjf^jf) + rn 2 ^. 
We use these formulae to calculate m and &HM for the various models discussed above and give 
the results in Table 2. 

2.4 Standard CDM and Variations 

We calculate Sh, m, and cfhai for a variety of models specified by their power spectra P(k). We 
begin by considering a family of models where the initial power spectrum is an n = 1, Harrison- 
Zel'dovich power law: Pi(k) = Ak. The power spectrum today is given by P, times an appropriate 
transfer function: P(k) = AkT 2 (k). Following Efstathiou, Bond, & White (1992) we choose the 
following parametric form for T(k): 



where a = 6Ah~ Mpc, b = 3.0/i Mpc, c = 1.7h~ Mpc, and v = 1.13. The power spectrum 
is therefore specified by T which determines the shape and A, or equivalently erg which sets the 
normalization. 5h oc erg and so, in Figure 1, we show 5h as a function of T /a%. The results are 
given for h = 0.5 and three different values of R. As expected, 5h increases with decreasing T. 

The Lauer & Postman (1994) velocity data on Abell clusters shows an large scale bulk motion 
that may indicate power on very large scales in excess of what is expected in most COBE normalized 
models. While there is still some controversy over the interpretation of this data (Strauss et al 
1994, Feldman & Watkins 1994) it is useful to consider the implications these results would have 
for determinations of the Hubble constant, the idea being that a universe with large scale bulk 
flows might also have large breathing mode fluctuations in the peculiar velocity field that foil our 
attempts to measure H. Jaffe & Kaiser (1995) apply a likelihood analysis to the Lauer & Postman 
(1994) data using Eq (|0]), and find a peak in the likelihood function at a$ ~ 0.3 and T ~ 0.025. 




(15) 






(16) 
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From Figure 1, we see that this implies a rather modest Sh- $h(R = 100 Mpc) = 0.05. They 
note however, that there is a rather wide range of parameters that are consistent with the Lauer & 
Postman data. We see from their Figure 1, for example, that V = 0.35 and a% = 1.2 is within the 
68% confidence intervals for the COBE quadrupole measurement and also within the 50% likelihood 
contour for the Lauer & Postman data. For these parameters, Sh — 0.12. 

We next consider the following models: (1) Standard CDM (sCDM) with COBE normalization 
(Bunn, Scott, & White 1994): This is essentially the model considered in Section 2.2 but the 
as = 1.38. (2) Mixed dark matter (MDM): We use Holtzman's (1989) linear power spectrum for a 
70% cold, 30% hot mixed dark matter model. Normalization is the same as in Model 1. (3) Tilted 
primordial P(k) with spectral index n > 1. Models with both CDM and MDM are used. Again we 
normalize to the COBE results on large scales following Bunn, Scott, & White (1994). The results 
for these models are summarized in Table 2. 









Tilted 


CDM 


MDM 




sCDM 


MDM 


n=1.15 


n=1.5 


n=1.5 


<$h(100 Mpc) 


0.074 


0.072 


0.098 


0.17 


0.16 


5 H (150 Mpc) 


0.040 


0.041 


0.052 


0.083 


0.083 


5h(200 Mpc) 


0.025 


0.026 


0.032 


0.049 


0.050 


(SH/H) vs. {5M/M) 












m(150 Mpc) 


-0.60 


-0.62 


-0.60 


-0.57 


-0.59 


ct(150 Mpc) 


0.10 


0.10 


0.13 


0.18 


0.10 



Table 2: Linear-theory 5h(R), m, and a for different power spectra P(k). 

Evidently, the rms peculiar Hubble flow on 100 Mpc scales can be quite substantial, especially 
for n > 1. These so-called blue primordial spectra (Lucchin et al. 1995) were introduced to boost 
power on 100 Mpc scales relative to COBE scales. Indeed, it is worth noting that the constraints 
on n from the COBE DMR experiment are not very restrictive, and n = 1.5 is easily allowed by 
the data (Bennett et al. 1994). A CDM model with n > 1 will no doubt lead to excessive power on 
small (< IMpc) scales and this has motivated Lucchin et al. (1995) to consider n > 1 with MDM. 
Another potential difficulty with these models is that they can lead to large temperature fluctua- 
tions in degree-scale CBR experiments. Consider the usual decomposition of the CBR anisotropy: 
ST/T(8, <p) = J2im a imYi m (0, 4>). For Gaussian perturbations, the temperature fluctuations are com- 
pletely specified by the variance C\ = (|a; m | 2 ). Following Kosowsky & Turner (1995) we estimate 
the effects of tilting the spectrum to be 

, ln(i/2)(n-l) 

C,(n)~- d(n = l) (17) 
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For I = 100 we find that the temperature fluctuations (oc C, ) are increased by a factor of 3 
over the n = 1 model for n = 1.15 and 50 for n = 1.5. The latter case is clearly unacceptable 
unless degree-scale fluctuations are damped by a period of early reionization. In any case, the next 
generation of CBR experiments will lead to tighter constraints on the primordial spectrum and 
help settle this issue. 

3 CBR Fluctuations from Voids 

We have seen that the possibility of living in a large underdense region opens up the prospect 
for solving the Ho/to problem without giving up spatial flatness or resorting to a cosmological 
constant (Suto et al. 1994; Wu et al. 1995; Moffatt & Tatarski 1994). In the simplest scenario, 
we would be at the center of an underdense region containing all of the objects used to measure 
H. We could then neglect safely the dynamics of the wall that separates this region from the rest 
of the Universe. For the CFHT and HST results this region would include the Coma cluster as 
well as clusters used to determine the peculiar velocity of Coma (as in, for example, Han & Mould 
(1992)). The region would therefore have to have a radius significantly larger than 100 Mpc. An 
even larger region would be required if one is to accommodate results from type II supernovae (e.g., 
Schmidt et al. (1994) which draw from objects up to 180 Mpc away. On the other hand, it may be 
possible to accommodate the CFHT and HST results with a local underdense region whose radius is 
comparable to 100 Mpc. This could arise if the Local Group and/or clusters such as Coma could lie 
on the surface of a void or perhaps on a ridge separating merging voids. Using N-body simulations 
to study interacting voids, Dubinski et al. (1993) find significant peculiar velocities tangential to 
the walls that separate merging voids, peculiar velocities that might help explain an anomalously 
high measurement of H. However, this scenario is inherently more complicated that the scenario 
in which the void contains the clusters used to determine H. In particular, the large scale peculiar 
velocity field is no longer characterized by a simple breathing mode. 

The mean density in our region of the Universe must be fairly small if local measurements of 
H are to be significantly higher than the global value. In Figure 2 we show Hl, the expansion rate 
inside a uniform region, as a function of the density in that region, f2 V oid- For |f2 V oid — 1| <C 1, we 
have the linear result discussed above: Hl — Hq = ^Hq (1 — S7 vo id)- The HST and CFHT results 
with quoted error bars are also shown. We see that to > 12Gyr and Hi > 63 km sec -1 Mpc -1 (the 
la lower bound for the HST result) imply a local density Slvoid 0.55, i.e., a nonlinear void. 

Voids somewhat smaller than the ones required here appear to be quite common. Redshift 
surveys reveal a network of voids with typical diameters of 50 — 60 Mpc (Kirshner et al. 1981; de 
Lapparent et al. 1986; Geller & Huchra 1989; da Costa 1991). There is also evidence for larger 
voids (Einasto, Joeveen, & Saar 1980, Bahcall & Soneira 1982) though, as one might expect, the 
larger voids are not nearly so empty. Of course, redshift surveys give the distribution of visible 
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galaxies and not the underlying matter distribution and so it is difficult to say what the true density 
contrast is inside the voids. Nevertheless, a picture emerges of a hierarchy of voids in which the 
density contrast decreases with increasing scale. 

Galaxy counts and redshift surveys offer the possibility of directly probing the density and 
velocity fields in our region of the Universe. However, they are each plagued by a number of 
difficulties. In particular, evolutionary effects can mask or mimic density inhomogeneities in the 
galaxy counts. In addition, there is the question of how one normalizes this data. On the other 
hand, attempts to map the peculiar velocity field must contend with sparse sampling; the rich 
clusters used to probe the velocity field are just too rare. 

Observations of anisotropies in the CBR present an alternative, and potentially clean probe 
of large scale density inhomogeneities. In this section we focus on constraints of large scale voids 
from CBR experiments. Previous attempts along these lines have used a statistical approach 
(Blumenthal et al. 1992; Piran et al. 1993). The idea was to construct a power spectrum that 
leads to a "reasonable number" of voids with the desired size and density contrast and to then see if 
this power spectrum is consistent with CBR observations. The discussion at the end of this section 
follows this approach. Here we focus on temperature fluctuations from individual voids and the 
primordial fluctuations that gave rise to them. In particular, we consider three distinct sources of 
temperature fluctuations: (1) primeval density fluctuations at the last scattering surface (primeval 
voids); (2) evolving voids between us and the last scattering surface (intermediate voids); and (3) 
the void enclosing us. 

3.1 Formalism 

Consider a spherically symmetric density perturbation in a matter dominated universe. The metric 
can be written (e.g., Kramer et al. 1980) 

ds 2 = -dt 2 + e x ^d X 2 + r 2 ( X , t)dfl 2 . (18) 

It is convenient to define the (spatially dependent) scale factor a(x,i): 

r( X ,t) =a( X ,t)x ■ (19) 

Let aj(x) = a(x,ti) and be the scale factor and background density at some initial time U and 
let 5(x) be the initial average density contrast enclosed by the x = constant surface: 

S(X) = Pi(x)/Pbi - 1- (20) 

where ~p~i(x) is the average density contrast enclosed by the x = constant surface. The equation of 
motion for a(x, t) follows from the Einstein equations: 

a(x,t) _ 1 +J(x) df(n) 
Oi(x) 25(x) dry 
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where /(r/) = rj — sin[h]ry for 5 > [<]0 and Hi(x) = (a/a)]*^. 
For a photon passing through the origin 



dx_ ± Vl-(4/9Wx)xV*g 
dt r' 

(Fang and Wu 1993) where the plus (minus) sign is for a photon moving away from (towards) the 
origin. The frequency shift of the photon is 

« = exp ( f S hdt) = exp ( [ tf *£±idt) , (23) 



Vf \Jti 2 / \J U xa' + a 

where Vi is the frequency of the photon emitted at t{ and Vf is the frequency of the photon observed 
at tt. Here and throughout, we use over dot to denote partial derivative with respect to t and prime 
to denote partial derivative with respect to x- The resultant temperature fluctuation for a blackbody 
distribution of photons passing through the center of the density fluctuation, in a spatially flat, 
matter-dominated universe, is 

rrp ..flue. bkgd / \ 2/3 , ■/ . - x 

By our convention, a negative 5T/T implies a cold spot in CBR sky. 

The frequency shift of a photon passing through the center of a void is found by numerically 



integrating Eq (24). The calculation is straightforward when the density perturbation is linear but 
requires a special approximation, to be discussed in section 3.3, once the void becomes nonlinear. 
Spherical symmetry is assumed for the voids through out our calculation, unless explicitly stated 
otherwise. 

3.2 Primeval Voids at the Last Scattering Surface 

We assume that the voids present today result from the gravitational amplification of primordial 
perturbations. The alternative would be to have the voids produced after recombination by some 
non-gravitational process. The evolution of the spacetime inside a void is calculated from Eq (pi]). 
For definiteness, we assume a density profile of the form: 

( #0, X < Xi- ) 

S(x) = { Soxl/ixl - X 3 i), Xi < X < Xo- \ (25) 
I 0, X >Xo- J 

so that the ratio xi '■ Xo determines whether one has a thick or thin wall. In Figure 3, we plot 5q 

as a function of f2 vo id today. We see, for example, that O voi( j < 0.55 requires 5o < —8.5 x 10~ 4 at 

recombination. 



11 



The correspondence between the size of the primeval void and its size today is trickier. In 



the linear regime, the size of the void can be easily calculated from Eq (21). In the nonlinear 

regime, matter from both inside and outside the void piles up on its surface (Maeda & Sato 1983a; 

Berschinger 1985). In any case, the size of the void can grow no faster than the expansion rate inside 

the void: A void that has £l vo id = 0.55 today grew no more than 20% (in comoving coordinates) 

from the time of recombination. 

Density perturbations at the time of recombination lead to temperature fluctuations in the CBR 

sky. If the scale of the perturbation I is much larger than the horizon size at recombination, ct Tec , 

the temperature fluctuations will be dominated by the Sachs- Wolfe effect (Sachs & Wolfe 1968). 

The expected temperature fluctuation is 

5T 5p ( I \ 2 
T p \ct iec J 

where we have omitted a pref actor of order 0.1. For I < ct rec , other effects such as the Doppler 
peak and adiabatic fluctuations in baryons become important (Hu, Sugiyama, & Silk 1995). 

Figure 4 shows our numerical calculation of the Sachs- Wolfe effect for a primeval void with 
5q = —8.5 x 10~ 4 and %l = 67/i _1 /(l + z rec ) Mpc and three different density profiles, xlss is the 
distance from the last scattering surface to the center of the void with negative values corresponding 
to the situation in which the last scattering surface lies behind the center of the void. In these 
calculations, we neglect the thickness of the last scattering surface since it is unimportant so long 
as Xi ^ 10/i _1 /(l + z rec ) Mpc. Note that the temperature fluctuations scale linearly with 5q. 

We see that sign and magnitude of the temperature fluctuation depends sensitively on where the 
last scattering surface intercepts the void. In addition, while the maximum temperature fluctuation 
is relatively insensitive to the thickness of the wall, the range in values of xlss ° ver which 5T/T is 
large increases as we increase the thickness of the wall. 

We wish to calculate the number of CBR fluctuations as a function of angular diameter and 
amplitude. Let n™£ d be the number density of voids that lead to CBR fluctuations with the desired 
amplitude and angular size. The expected number of fluctuations seen by a given experiment will 
be 



\ i- ~r z rec J 



14/i~ 3 Mpc 3 



Void 



Xi 



V (27) 



.100Mpc/(l + z T , 

The COBE DMR experiment has mapped the entire sky with 7° resolution. Primeval voids 
with xi ^ 350/i -1 /(l + z rec ) and <5o < —8.5 x 10 -4 lead to fluctuations that would have been easily 
seen in this experiment. We therefore require (Af) < 1 for these very large voids. This implies a 
tight constraint on their number density today, n° oid = i^ d /(l + z Tl v ' 



j rcc J 



"void ^ 2xlQ- 11 h 3 



100/(1 + z rcc )Mpc 



Xi 



Mpc -3 
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< 5 x 1CT 11 MpcT 3 . (28) 

Degree-scale measurements of temperature fluctuations in the CBR place constraints on smaller- 
scale voids. The MAX and MSAM experiments cover ~ 40 square degrees in total (V ~ 0.001) 
with 0.5° resolution and find an rms temperature fluctuation of \8T/T\ < 6 x 10~ 5 (Meinhold and 
Lubin 1991; Alsop et al. 1992; Meinhold et al. 1993; Gundersen et al. 1993; Devlin et al. 1994; 
Clapp et al. 1994; Cheng et al. 1994). Figure 5 shows the absolute value of 5T/T in the direction 
of the center of the void, averaged over —0.5 < Xlss/Xo < 0.5. Prom Figure 5, we see that primeval 
voids with Sq < —8.5 x 10~ 4 , xi k 50/i~ 1 (l + z rC c) 1 Mpc, and xo : Xi = 3 : 1 (thick wall case) 
or Xi k 80/i _1 (l + z rcc )~ l Mpc and xo : Xi = 9 : 8 (thin wall case) can easily lead to temperature 
fluctuations that exceed the observational limit. Several individual patches with \5T/T\ as large as 
1 x 10~ 4 have been observed though it is not as yet clear if these rare large fluctuations represent true 
CBR anisotropies or foreground contamination. In any case, primeval voids with 5o < —8.5 x 10~ 4 
and xi k, 50/i -1 (l + 2; rec ) _1 Mpc (thick wall) or xi k, 80/t _1 (l + z^)' 1 Mpc (thin wall) should be 
rare in the currently covered area of MAX and MSAM. To be conservative, we require (AT) < 10. 
This yields the constraints 

n° mid < 3 x 10- 7 /i 3 Mpc- 3 . (29) 

This implies that < O(0. 1) of the volume of our universe can be occupied by these voids. 

3.3 Voids at Intermediate Redshift 

Temperature fluctuations due to nonlinear structures between us and the last scattering surface 
were first considered by Rees and Sciama (1968). For subhorizon-sized structures, the temperature 
fluctuation will be 0((R/ct) 3 ) where R is the physical size of the object at the time t when the 
photon crosses it. 

We use the formalism of Section 3.1 to calculate temperature fluctuations due to an intermediate 
nonlinear void of arbitrary interior density. For simplicity, we assume the a thin wall separates the 
void from the rest of the Universe. For the special case of a vacuum void, Thompson and Vishniac 
(1987), and Martinez-Gonzalez et al. (1990), find 

5T fR\ 3 T /8 16 16 2 \ 

: cos * -7 - — - — cos 2 f (30) 



T \ctj V9 27 81 

where \E f is the angle between the photon's direction and a line from the center of the void to the 
exit point of the photon. 7 describes the expansion of the void: Rat 1 . 

The numerical integration of Eq ( |24| ) must be modified in order to handle the nonlinear voids 
considered here. In the formalism outlined above, the initial density perturbation is divided into 
concentric, spherical shells. However, once the density fluctuation becomes nonlinear and matter 
starts to pile up on the wall of the void, shell-crossing will occur. To avoid complications due 
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to shell-crossing, we choose instead to divide the density profile into concentric shells at the time 
when the photons exit the void. The "evolution" of the void is now found by solving the same 
equations as before, but letting time run backwards. However now, there is no shell crossing and 
the frequency shift of the photons can be calculated using the same methods as in section 3.1. 

The expansion rate inside the void is still calculated by evolving the initial perturbation (i.e., 
the perturbation at recombination) forward in time. This calculation not only determines Hl but 
also the velocity perturbation associated with the void at the time when the photons exit the void. 
The latter is required for our backward-time integration. The index 7 for the expansion of the 
wall is put in by hand. In general, 0.667 < 7 < 0.8. 7 = 2/3 corresponds to the case where the 
wall expands with the background, a situation that will arise when the wall collides with the walls 
enclosing neighboring voids. 7 = 0.8 corresponds to a single isolated void. The result is easily 
derived from either energy or momentum conservation (Maeda and Sato 1983a; Berschinger 1985). 

For a nearby vacuum void with a thin wall and a radius of 100/i _1 Mpc, eq. ( p0| ) yields 
approximately —1 x 10~ 5 for 7 = 0.8, and —2.5 x 10~ 5 for 7 = 2/3 in the direction of the center 
of the void. Our calculations yield —7 x 10 -6 and —2.6 x 10 -5 respectively. Our calculations also 
confirmed the scaling law of 5T/T oc (R/ct) 3 . Empirical comparisons indicate that our results are 
reliable to order 0(1O -6 ). 

Figure 6 shows our results for a void with a lOO/i -1 Mpc radius and with different ^ V oid under 
a thin wall approximation. 5T/T is the CBR temperature fluctuation in the direction through the 
center of the void. It is assumed that we are located just outside the wall. The figure indicates 
that the ST/T due to a nearby 100/i _1 Mpc void is not large enough to be constrained by COBE 
measurements. However, according to the (R/ct) s scaling law, voids with a radius > 200/i~ 1 Mpc 
(thus yielding > 8 times more temperature fluctuation) and a sufficiently low O v oid can be in conflict 
with the COBE observation. 

3.4 The Void Enclosing Us 

For spherical voids between us and the last scattering surface, the leading order (O(Rfdo)) and 
the second order (0(R/ 'do) 2 ) contributions to the CBR anisotropy cancel (Rees & Sciama 1968). 
For the void enclosing us, this cancellation does not occur and we have a contribution to the dipole 
component of the CBR of 0(R/do)- This contribution can be viewed simply as arising from our 
peculiar motion relative to the CBR rest frame and its amplitude is easy to estimate: 



Here r is the our distance to the center of the void, as illustrated in the geometry of Figure 7. 
A detailed calculation, making use of the formalism developed above, leads to similar results. In 
Figure 8, we show the amplitude of the dipole moment as a function of r for < S^ vo id < 0.55. 




(31) 
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The requirement that \5T/T\ < 10 3 leads to the constraints r < 7h x Mpc for a vacuum void, and 
r < 22h- 1 Mpc for ft void = 0.55. 

The contribution to higher order moments from the void enclosing us is 0((R/ ct^f) and can 
lead to unacceptably large fluctuations for R > 200/i -1 Mpc, unless we are essentially at the center. 

Thus far we have assumed perfect spherical symmetry for the voids. This is no doubt unrealistic, 
though it does appear that voids become more spherical with time (Blaes et al. 1990). Clearly, 
deviations from spherical symmetry will lead to CBR distortions. Consider, for example, a bump in 
the wall with height AR (Figure 7). We estimate of its order of magnitude as follows. The deviation 
from spherical symmetry should contribute to temperature fluctuation at orders {R/ cto) m (AR/ R) n , 
where m and n are positive integers. The leading order, with m = n = 1, is simply a dipole 
resembling eq. (|3l|). The next order terms, with m + n = 3, contribute to higher order moments. 
(For spherical symmetry (Rees and Sciama 1968), only the m = 3, n = term is relevant.) By 
requiring contributions at these orders to be smaller than 10 -5 , the level of anisotropy measured 
by COBE at > 10° angular scale, we can constrain the size of the deviation AR. For R < 100/i _1 
Mpc, a bound is estimated by considering the m = 1, n = 2 term: 



5T 



T 



R \ /AR 



(A "M^J (it) * 10 " (32) 

where AH is the difference between the expansion rates inside and outside the void, which is 
roughly Hq(1 — r2 V oid)/3 in linear theory. Eq. (|32"|) amounts to 



2 



AR ^ / lOO/i^Mpc , . 

at > 10° angular scale. If R > 100/i -1 Mpc, a bound, estimated from the m = 2, n = 1 term, gives 



AR 0.02 / 100/i _1 Mpc\ 



2 



R ~ 1 — il VO id \ R ) ^ ^ 

Therefore, a 200/i _1 Mpc radius void must be spherically symmetric to ~ 1% to avoid a contradic- 
tion with the COBE results. 



4 Discussion and Conclusions 

The Hubble constant sets both the distance and time scales for the Universe and therefore plays 
a central role in all cosmological models. Numerous observations place its value between 65 and 
95kmsec~ 1 Mpc -1 . If indeed the true value is in this range, and if the Universe is as old as the 
globular cluster experts say it is, then fairly radical modifications to the standard f2 m = 1 paradigm 
will be called for. 

In this paper, we explore the possibility that the true Hubble constant Hq ~ 50 km sec -1 Mpc -1 , 
not because of systematic errors in current observations, but because these observations do not reach 
deep enough into the cosmos. 



15 



Possible discrepancies between the locally measured value for H and the true value are quantified 
in 5h = 5h(R) where R is the radius of the sample volume. We outline a simple procedure for 
calculating 5h based on linear perturbation theory. Where applicable, we compare the results with 
those based on the N-body simulations of TCO. Our conclusions are as follows: 

• Linear and nonlinear theories are in good agreement on large scales. 

• For R > 100 Mpc, 5h is insensitive to the type of dark matter present in the model. 

• 5h is correlated with the mass fluctuation within the sample. In particular, a measurement 
of 5M/M leads to the prediction 5H/H = -(0.6 ± 0.1S)SM/M. 

• Sh(R = 100 Mpc) ~ 0.07 for standard n = 1 CDM as well as n = 1 MDM. If these models 
describe our Universe, then the HST and CFHT results, in all likelihood, reflect the true value. 

• 5h(R = 100) Mpc > 0.1 for nonstandard models with extra power on large scales. These 
include an n = 1 CDM-like model with T = 0.35 and a% = 1.2 (see Eq ([TB])) and either CDM or 
MDM with a tilted, n = 1.5 primordial spectrum. If our Universe is described by either of these 
models, then measurements of H will have to reach ~ 200 Mpc where 5h drops well below 0.1. 

Galaxy counts and redshift surveys allow us the opportunity to directly probe the density and 
velocity field in our region of the Universe. Unfortunately, the interpretation of these observations 
is not so straightforward. In the case of the galaxy counts, for example, it is difficult to separate 
density inhomogeneities from evolutionary effects. On the other hand, analysis of the velocity 
data is made difficult by sparse sampling. CBR anisotropy measurements offer an alternative, and 
potentially clean probe of the large scale inhomogeneities discussed in this paper. Here, we calculate 
the expected temperature fluctuations due to voids at the last scattering surface; between us and 
the last scattering surface; and by the void enclosing us. The results are compared with existing 
data from the COBE DMR experiment, and by degree-scale experiments such as MAX and MSAM. 
A summary of our conclusions follows: 

• No more than a handful of voids with radii > 400h _1 Mpc can exist in our Hubble volume. 

• Smaller voids are constrained by degree-scale CBR experiments. We find that the number 
density today of voids with radii > 60h -1 Mpc and < 0.55, under the current sky coverage of 
degree-scale CBR experiments, is < 3 x 10~ 7 /i 3 Mpc -3 , occupying O(0. 1) of our Hubble volume. 

• If the void enclosing us has a radius > 200h -1 Mpc, large dipole and quadrupole anisotropies 
will be induced, unless, of course, the void is very nearly spherical, and we are near the center. 
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Figure Captions: 



Figure 1. Sjj, the rms deviation in the Hubble flow, as a function of the transfer function param- 
eterized by T (Eq (|l6|)) for scales R = 100, 150, and 200 Mpc, h = 0.5, ft m = 1, and a 8 = 1. 5 H 
scales linearly with og. 

Figure 2. The relation between Hl, £l yo id, and to m an = 1, A = universe. 

Figure 3. The relation between the initial density fluctuation 5q in a primeval void and the density 
parameter £l vo id it would have today. 

Figure 4. The temperature fluctuation of CBR in the direction of the center of a primeval void, 
due to the Sachs- Wolfe effect, vs. the position of the last scattering surface relative to the center of 
the void. The initial density fluctuation is assumed to be 5q = —8.5 x 10 -4 . Different lines stands 
for cases with different wall thickness. 

Figure 5. The averaged temperature fluctuation magnitude in CBR from a primeval void, vs. 
the coordinate radius of the void, for three different wall thickness. The dot-dash line shows the 
temperature fluctuation vs. the coordinate radius of the void when the last scattering surface cut 
through the center of the void. The initial density fluctuation is assumed to be 5q = —8.5 x 10~ 4 . 

Figure 6. The temperature fluctuation magnitude in CBR from a intermediate void with a radius 
of lOO/i -1 Mpc, vs. the density parameter f2 VO id inside the void. We are assumed to be just outside 
the void. The solid line: 7 = 2/3; the dash line: 7 as calculated from momentum conservation of 
the wall. 

Figure 7. The geometry of a void enclosing us. The dashed line illustrate a deviation of the wall of 
the void from a perfect sphere. AR is the height of the deviation. 

Figure 8. The amplitude of a CBR dipole moment \5T /T\&\ vo \ e contributed by a void enclosing us 
vs. the distance between us and the center of the void r. No peculiar velocity other than that due 
to the faster expansion inside the void is assumed for us. 
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